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We have measured with a tunnel probe the energy distribution function of 
quasiparticles in silver diffusive wires connected to two large pads ("reser- 
voirs"), between which a bias voltage was applied. From the dependence 
in energy and bias voltage of the distribution function we have inferred the 
energy exchange rate between quasiparticles. In contrast with previously ob- 
tained results on copper and gold wires, these data on silver wires can be well 
interpreted with the theory of diffusive conductors either solely, or associated 
with another mechanism, possibly the coupling to two-level systems. 
PA CS numbers: 73.23.-b, 73.50.-h, 71.10. Ay, 72.70. +m. 



The present understanding of metals at low temperature relies on Lan- 
dau's theory of Fermi liquids. In this theory, the elementary excitations 
of the electron fluid are nearly independent fermionic quasiparticles.El The 
residual interactions depend on the efficiency of the screening of Coulomb 
interactions, and increase if electrons are scattered by impurities, surface 
or lattice defects.B'0 These interactions can be probed through the shape 
of the energy distribution function in an out-of-equilibrium situation. We 
have found in previous experiments on coppeiB and goldi wires that the 
energy exchange rate between quasiparticles is stronger and has a different 
energy dependence than predicted by Altshuler and Aronov's (AA) theory 
of diffusive conductors.! In the present paper, we report measurements of 
distribution functions in silver wires, where the interactions appear to be in 
closer agreement with AA theory. 

The experimental setup used to probe the energy exchange rate is shown 
in Fig. pi. A metallic diffusive wire of length L is connected at both ends 
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Fig. 1. Experimental layout: a metallic wire of length L is connected to large 
reservoir electrodes, biased at potentials and U. In absence of interaction, 
the distribution function at a distance X = xL from the grounded electrode 
has an intermediate step / (E) = 1 — x for energies between — eU and 
(solid curves) (we assume U > 0). When interactions are strong enough to 
thermalize electrons, the distribution function is a Fermi function, with a 
space-dependent temperature and electrochemical potential (dotted curves). 
In the experiment, the distribution function is obtained from the differential 
conductance dI/dV(V) of the tunnel junction formed by the wire and a 
superconducting electrode placed underneath. 



to large and thick electrodes called "reservoirs" in the following. The quasi- 
particle energy distribution f(x,E) at a distance X = xL from the right 
electrode, is obtained from the differential conductance dl/dV (V) of a tun- 
nel junction between the wire and a superconducting electrode.oi A voltage 
difference U is applied between the reservoirs in order to implement a sta- 
tionary out-of-equilibrium situation. The shape of the distribution function 
f(x, E) depends on the average number of inelastic collisions a quasiparti- 
cle experiences during its diffusive motion from one of the electrodes to the 
position X, and on the amount of energy exchanged at each collision. 

We present the results of four experiments, labeled A, B, C and D in 
the following. The wire length L = 5, 10, or 20 /im is appended to the label; 
for instance the samples A5 and C20 are 5 and 20 /im-long, respectively. All 
samples were fabricated by electron-gun evaporation of silver at several an- 
gles through a PMMA suspended mask patterned using e-beam lithography. 
The substrate was, as in our experiments on Cu and Au, thermally oxidized 
silicon. Samples with same labels (B5 and BIO on the one hand, D20a and 
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Fig. 2. Measured distribution functions for U = 0, 0.1 and 0.2 mV in 
samples B5, BIO, D20a and D20b (see Table 1). In the top left panel, the 
dotted line is the prediction for the non- interacting regime (Eq. (2)) for 
U = 0.2 mV. 



D20b on the other hand) were fabricated simultaneously. The thickness of 
the wires is 45 nm and their width w ranges between 65 nm and 150 nm. The 
electrodes at the ends of the wires are 400 nm-thick silver pads with an area 
of about 1 mm 2 , thereby implementing adequate reservoirs J3 The supercon- 
ducting probes, made of aluminum,! were positioned at x = 0.5 and/or at 
x = 0.2. The areas of the tunnel junctions are w x 150 nm and their tunnel 
resistances range between 23 kfi and 130 kfi. We estimated the diffusion co- 
efficient D of quasiparticles, and hence the diffusion time td = L 2 /D, from 
the low-temperature resistance of the wires. The samples were mounted in 
a copper box thermally anchored to the mixing chamber of a dilution refrig- 
erator. Electrical connections were made through filtered coaxial lines p and 
measurements were carried out at a temperature of 40 mK. 

The distribution functions obtained at U = 0, 0.1 and 0.2 mV are shown 
in Fig. |2[ For [7 = 0, f(E) is close to the expected Fermi function at the 
temperature of the refrigerator. For U ^ 0, the functions f(E) measured 
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near the middle of the wire B5 (top left panel) display a sharp double step 
with a plateau at height 1 — x = 0.57. For comparison, we have plotted as a 
dotted line the best fit at U = 0.2 mV with a linear combination of the Fermi 
functions of the reservoirs ,0 which is the expected distribution function in 
the non-interacting regime (see Refs.B0 and Eq. (||) below). The deviations 
from this regime are more apparent in sample D20a (bottom left panel), for 
which the diffusion time td = 20 ns is significantly longer than in sample 
B5 where td = 1-8 ns. In the right panels of Fig. |2| we show the distribution 
functions measured at the lateral position of samples BIO (top right panel) 
and D20b (bottom right panel). The height of the plateau well agrees with 
1 — x ~ 0.8, and the distribution functions are again more rounded for the 
wire with the longer diffusion time. 

The distribution functions obtained in the experiments on copper and 
gold wires systematically display a scaling propertyja as shown for a copper 
sampled in the right panel of Fig. ||, f(x, E) only depends, at each position, 
on the reduced variable E/eU. Such a scaling law is not observed in our silver 
samples: in particular, the slope of the plateau in the distribution functions 
of the wire D20a increases with U when plotted in reduced units (see the 
left panel of Fig. ||). This indicates that interactions between quasiparticles 
have a different energy dependence in the two types of samples. Interactions 
are also weaker in silver samples: as illustrated in Fig. [3|, similar distribution 
functions are obtained for longer diffusion time in silver than in copper or 
gold samples. 

In order to compare the energy distribution functions f (x,E) obtained 
experimentally with the theoretical predictions, we now explain how the 
energy exchange rate between quasiparticles determines / (x, E) . The sta- 
tionary distribution function / (x, E) obeys the Boltzmann equation: 

1 °' J U ' L:) +l co ii(x,E,{f}) = (1) 



td dx 2 



where I co n (x, E, {/}) is the collision integral due to interactions between 
quasiparticles. The boundary conditions are imposed by the reservoirs at 

both ends: / (0, E) = (l + exp ^) ^ and / (1, E) = (l + exp ^f) ~* . 

In the absence of inelastic scattering T co \\ = and the distribution 
function f (x, E) is:0 

f (x,E) = (l-x)f(0,E)+xf(l,E). (2) 

The function /q (x, E) has a well-defined plateau for \eU\ 3> ksT, as observed 
in sample B5 (see Fig. 2). 

The collision term I co \\ (x, E, {/}) is the difference of two terms: an in- 
collision term, the rate at which particles are scattered into a state of energy 
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Fig. 3. Distribution functions for U = 0.1, 0.2, 0.3, and 0.4 mV, plotted as 
a function of the reduced energy E/eU. Left panel: Ag sample D20a; right 
panel: Cu sample, L = 5 /im. 

E, and an out-collision term: 

Zeoii (x, E, {/}) = Z£u (*, {/}) - Z£S (x, tf, {/}) (3) 

with 

l^(x,E,{f}) = JdedE'K(e) (4) 

x /1+£,b(1 ~~ /l,E-e)/i'(l ~~ f%'+e) 

where the shorthand stands for / (x, . Following Landau's approach,^ 
we have first assumed that the dominant process is a two-quasiparticle in- 
teraction. Moreover the interaction is assumed to be local on the scale of 
variations of the distribution function. The kernel function K (e) is pro- 
portional to the averaged squared interaction between two quasiparticles 
exchanging an energy e. The scaling property observed for copper samples 
implies^ K (e) oc e~ 2 , as opposed to the AA prediction K (e) oc e -3 / 2 for 
a diffusive conductor in the ID regime.@ In silver samples we have assumed 
that the interaction kernel still obeys a power law K (s) = K a e~ a , with k& 
and q taken as fitting parameters. These best fit theory curves, obtained 
with the parameters given in the colums "fit K a e~ a " of Table 1, are plotted 
in Fig. |4| with full squares. For comparison, the best fits obtained with 
the exponent set at its predicted value a = 3/2 are plotted with open dia- 
monds. Note that the same fitting parameters were used for all the samples 
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Table 1. Samples parameters and fit parameters. The wire width w is in 
nm, the diffusion constant D in m 2 s _1 . The fit parameters n a are given in 
ns _1 meV a_2 , and ktls is given in ns _1 meV _1 . 



of a given experiment. We obtained the exponent values a = 1.6 ± 0.2 and 
1.4 ± 0.2 for experiments A and B, respectively, which are compatible with 
the prediction a = 3/2 (the diamonds in the top panels of Fig. |] are prac- 
tically superimposed with the squares). This compatibility is not found for 
the experiments C and D, for which a = 1.2 ± 0.15 : the quality of the fits 
is visibly degraded by imposing a = 3/2, as shown in the bottom panels of 
Fig. |4j. The slope of the intermediate plateau in the theoretical curve (open 
diamonds) is systematically too large at U = 100 fiV and U = 200 /uV. 

In order to explain the discrepancy between the predicted exponent 
a = 3/2 and the value a = 1.2 ± 0.15 extracted from experiments C and 
D, we have investigated theoretically the effect of two-level systems (TLS) 
on the energy exchange between quasiparticles. The relevance of_TLS on 
phase relaxation has recently been suggested by several authors .EJ'E] We 
assume here that the quasiparticles are weakly coupled to the TLS, which 
are equally distributed along the wire. These TLS could be atoms moving 
across crystalline defects, or impurities in the crystal, for example. The 
TLS are assumed to have a flat energy distribution, and to be all equally 
coupled to the quasiparticles. We treat the absorption and emission of en- 
ergy by the two-level systems in the perturbative limit: the absorption rate 
r^_(e) (respectively the emission rate r^(£:)) by a TLS at position x with an 
energy separation e is given by Fermi's golden rule: r5.(e) = Xp^(e)h x (£) 
(resp. r^(e) = \p+(e)h x (— e)). In these expressions, A is the coupling 
constant between quasiparticles and the TLS, p x _(e) (resp. p+(e)) is the oc- 
cupation probability of the low-energy level (resp. the high-energy level), 
and h x (e) = J dE f%(l — f%- e )- Assuming that the TLS reach a lo- 
cal equilibrium with the quasiparticles, i.e. r+(e) = r^(e), one obtains 
(e) = h x {±e)/{h x {-e) + h x {e)). Each term Ap_(e)/ B (l - f E -e) in Y + (e) 
corresponds to an energy transfer of e from a quasiparticle at energy E to a 
TLS. Therefore, it gives an "out" collision term in the Boltzmann equation 




Fig. 4. Continuous lines in all four panels: measured distribution functions, 
plotted as a function of the reduced variable E/eU for U = 0.1, 0.2, 0.3 and 
0.4 mV, for the same samples as in Fig. 2. Successive curves were shifted 
vertically by 0.2, for clarity. Full squares are the best fits of the data to 
the solution of the Boltzmann equation with an interaction kernel K (e) = 
K a e~ a . Open diamonds are the best fits obtained with a set to the theoretical 
value 3/2, and adjusting the prefactor k 3 / 2 - Detailed fit parameters are given 
in Table 1. 
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(Eq. ((T|)) for quasiparticles at energy E. One then finds directly that the cou- 
pling to TLS can be written as an effective kernel function in Eq. (Q), which 
depends on the local distribution function: K x s (e) = KThs/(h x (—e) + h x (e)). 
The parameter ktls is proportional to the density of TLS and to the cou- 
pling constant A. The effective kernel K x s (e) is not a power-law function 
of e: at energies large compared to ell, K x s (e) oc e _1 , but K x s (0) remains 
finite. Without assuming an unrealistic heating of the reservoirs, it was not 
possible to fit the data with this model alone. However, we found that all 
the experimental data could be well accounted for by assuming the presence 
of two phenomena: direct quasiparticle-quasiparticle interaction, described 
by the AA theory:! K(e) = k 3 / 2 £ -3//2 ; and quasiparticle-TLS coupling, de- 
scribed by ifg ff (e). The parameters ktls and K3/2 for the best fits (which 
are not shown because they are hardly distinguishable from the fits with the 
best value of a in Fig. 4) are given in the columns "fits TLS+AA" of Ta- 
ble 1. Note that the error bars on ktls an d K3/2 are correlated: the weighted 
sum of ktls an d K3/2 must remain constant. Apart from sample A5, the 
theoretical value for are of the same order of magnitude as ^3/2- For 
the experiments A and B, extra contributions to AA theory are minimal. In 
contrast, in the experiments C and D the slope of the plateau in the distri- 
bution function is well explained by the TLS, whereas the AA mechanism, 
which dominates at low energy, is responsible for the rounding of the steps. 

The differences that we observed between copper, gold and silver wires 
in the energy exchange rate experiments were also found in measurements 
of the phase coherence time in wires fabricated using the same procedure. 
In a silver sample, the temperature dependence of follows closely the 
theoretical prediction down to the base temperature T = 50 mK of the 
refrigerator (where 70 = 9 ns), whereas saturates below 1 K at = 1 ns 
in our copper sample, and below 6 K at = 10 ps inpur gold sample. 
These measurements are further discussed in this volumetZI . 

In conclusion, we have found that interactions between quasiparticles in 
silver wires are much weaker and have a different qualitative behavior than 
in copper or gold wires. The energy exchange rate in our silver samples is 
close to the theoretical predictions for a diffusive medium in the ID regime, 
provided that one includes an extra contribution, which might be due to two- 
level systems. The difference in the behavior of interactions in gold, copper 
and silver samples, as seen from the energy exchange rate experiments, is 
correlated to the measured saturation of ta at low temperature which did 
occur in copper and gold, but not in silver. A possible interpretation could 
be that two- level systems are more numerous and/or better coupled to quasi- 
particles in copper and gold than in silver, leading to faster phase and energy 
relaxation. Moreover, preliminary calculations in the strong coupling regime 
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account for the scaling property of the distribution functions found in copper 
and gold samples.0 Experiments are in progress to test this interpretation. 
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